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Abstract 

We study the following charmless baryonic three-body decays of B mesons: B + — > pA7, B + — > 
pAtr and B° — > pAir~ . The partial branching fractions as a function of the baryon-antibaryon 
mass and the polar angle distributions of the proton in the baryon-antibaryon system are presented. 
This study includes the first observation of B + — > pA7r°, which is measured to have a branching 
fraction of (3.00lg;^ ± 0.33) x 10" 6 . We also set upper limits on branching fractions of the two- 
body decays B° -► pE* _ , B° -► A A, B + -► and B + -► A+A at the 90% confidence level. 
These results are obtained from a 414 fb -1 data sample collected near the T(45) resonance with 
the Belle detector at the KEKB asymmetric-energy e + e~ collider. 

PACS numbers: 13.40.Hq, 14.20.Dh, 14.40.Nd 
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After the first observation of charmless baryonic B meson decay, B + -> ppK + 0, 3, 
various three-body baryonic decays were found [l, 0, lEf. The dominant contributions for 
these decays are presumably via the b — > s penguin diagram as shown in Fig. [1] for the 
case of B + — > PA7. A common experimental feature of these decays is that the baryon- 
antibaryon mass spectra peak near threshold. This feature was conjectured in Ref. [f| and 
has recently aroused much theoretical interest 0]. Detailed information from the polar angle 
distributions [s[ and Dalitz plot 0] offer better understanding of the underlying dynamics. 

In this paper, we study the following three-body charmless baryonic decays of B mesons: 
B + — > pAj, B + — > pAir° and B° —>■ pA^~ . The partial branching fractions as a function of 
the baryon-antibaryon mass and the polar angle distributions of the proton in the baryon- 
antibaryon system are presented. It is interesting to compare the results with theoretical 



predictions lid, 1 1 if ] . Since the A hyperon could be a useful tool to probe the helicity selection 
rule for the b — > s process 12[ , we investigate the proton polar angular distribution from 
A decays. We also search for intermediate two-body decays in these three-body final states. 
This is motivated by the observations of two-body decays of charmed baryons [l3j]. Using 
topological quark diagrams for B decays and the assumption of SU(3) flavor symmetry, 
various two-b ody charmless baryonic decay modes should be observable with a data sample 
of -400 fb" 1 Q. 

We use a 414 fb -1 data sample consisting of 449 xl0 6 BB pairs collected with the Belle 



detector at the KEKB asymmetric energy e + e (3.5 on 8 GeV) collider [151 ] . The Belle 



detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter (ECL) composed of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return 
located outside the coil is instrumented to detect K\ mesons and to identify muons. The 
detector is described in detail elsewhere (l6| . The following two kinds of inner detector 
configurations were used. A 2.0 cm beam pipe and a 3-layer silicon vertex detector was 
used for the first sample of 152 xl0 6 BB pairs, while a 1.5 cm beam pipe, a 4-layer silicon 
detector and a small-cell inner drift chamber were used to record the remaining 297 xlO 6 BB 
pairs 
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The event selection criteria are based on the information obtained from the tracking sys- 
tem (SVD and CDC) and the particle identification system (CDC, ACC, TOF and ECL). 
They are optimized using Monte Carlo (MC) event samples produced by the EvtGen gener- 
ator [18] and GEANT [19 !i ] detector simulation. All primary charged tracks are required to 
satisfy track quality criteria based on the track impact parameters relative to the interaction 
point (IP). The deviations from the IP position are required to be within ±0.3 cm in the 
transverse (x-y) plane, and within ±3 cm in the z direction, where the +z axis is opposite 
to the positron beam direction. For each track, the likelihood values L p , L^, and that 
it is a proton, kaon, or pion, respectively, are determined from the information provided by 
the particle identification system. The track is identified as a proton if L p /(L P + L#) > 0.6 
and L p / (L p + L w ) > 0.6, or as a pion if L n / (L K + L w ) > 0.6. For charged particles with mo- 
menta around 2 GeV/c, the proton selection efficiency is about 84% (88% for p and 80% for 
p due to larger p cross sections) and the fake rate is about 10% for kaons and 3% for pions. 
Candidate A baryons are reconstructed from pairs of oppositely charged tracks — treated as 
a proton and negative pion — whose mass is consistent with the nominal A baryon mass, 
1.111 GeV/c 2 < M pn - < 1.121 GeV/c 2 . The A candidate should have a displaced vertex 
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FIG. 1: A b — * s penguin diagram for B + — > pA*y. 



and its momentum direction should be consistent with a A originating from the IP position. 
For particle identification of the A daughters (i.e. secondary charged tracks), we require 
only L p /(L p + L n ) > 0.6 for the proton, but do not impose any additional requirements on 
L p / (L p + Lk) for the proton or L n / [Lk + L n ) for the pion. Photon candidates are selected 
from neutral clusters in the ECL. Pairs of photons with invariant masses in the range 115 
MeV/c 2 < m 77 < 152 MeV/c 2 are used to form 7r° mesons. The measured energy of each 
photon in the laboratory frame is required to be greater than 50 MeV. The momentum of 
the 7T° in the laboratory frame should be greater than 200 MeV/c. The cosine of the decay 
angle should satisfy | cos 6^1 < 0.9, where 6 1 is the angle between the photon direction and 
the negative of the laboratory frame direction in the ir° rest frame. The primary photon 
from the B + — > pA'y decay must satisfy the following additional requirements: it should be 
in the barrel region (with polar angle between 33° and 128°) and have an energy greater 
than 500 MeV. We discard the primary photon candidate if, in combination with any other 
photon above 30 (200) MeV, its mass is within ±18 (±32) MeV/c 2 of the nominal mass of 
the 7T° (if) meson. 

Candidate B mesons are reconstructed in the B + — > pA7, B + — > pAir and B° — > pAir~ 
modes. We use two kinematic variables in the center-of-mass (CM) frame to identify the 



reconstructed B meson candidates: the beam energy constrained mass M^ c = y E^ eam — p 2 B , 
and the energy difference AE = EB — E^ eSum , where -Ebeam is the beam energy, and pb and Eb 
are the momentum and energy, respectively, of the reconstructed B meson. The candidate 
region is defined as 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.16 GeV < AE < 0.5 GeV 
for the 7r°/7 mode (—0.1 GeV < AE < 0.3 GeV for the 7r~ mode). The signal peaks in 
the subregion 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.135 GeV< AE < 0.074 GeV for 
the 7r°/7 mode (|A£J| < 0.03 GeV for the ir~ mode). The lower bound of AE is chosen to 
exclude possible contamination from so-called "cross-feed" baryonic B decays, i.e. four-body 
decays with a missed daughter. 

The background in the candidate region arises predominantly from the e + e~ — > qq (q — 
u, d, s, c) continuum. We suppress the jet-like continuum background relative to the 
more spherical BB signal using a Fisher discriminant [lo[ that combines seven event shape 



variables as described in Ref. 2JJ. The B + — > pAir mode has more background than the 
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other modes and therefore we add the missing mass to the Fisher variable. The missing 
mass is determined from the rest of the detected particles (treated as charged pions or 
photons) in the event assuming they are decay products of the other B meson. We form the 
signal (background) likelihood C s (£&) by combining probability density functions (PDFs) 
for the Fisher discriminant and the cosine of the angle between the B flight direction and the 
beam direction in the T(4S) rest frame. The signal PDFs are determined using signal MC 
simulation; the background PDFs are obtained from the side-band data with M^ c < 5.26 
GeV/c 2 . We require the likelihood ratio 71 = C S /(C S + Cb) to be greater than 0.75, 0.85, 
and 0.80 for the pA^, pA7r° and pAn~ modes, respectively. These selection criteria are 
determined by optimizing n s /^n s + rib, where n s and rib denote the expected numbers of 
signal and background events, respectively. We use the branching fractions from our previous 
measurements (a, [s| in the calculation of n s . The branching fraction of B + — > pAw® is 
assumed to be one half that for B° — > pAn~ [10]. If there are multiple B candidates in a 



single event, we select the one with the best TZ value. We apply a A+ — ► Att + veto for the 
B° — » pAir~ mode: candidate events with a reconstructed Atc + mass in the range 2.26-2.31 
GeV/c 2 are excluded. 

We perform an unbinned extended likelihood fit that maximizes the likelihood function, 



e -(N s +N b ) N_ 



J] ^N s P s (M hc ., AEi) + N b P h {M hCi ,AE i ) 



4=1 



to estimate the signal yield in the candidate region. Here P s (Pb) denotes the signal (back- 
ground) PDF, N is the number of events in the fit, i is the event index, and N s and Nb are 
fit parameters representing the number of signal and background events, respectively. 

For the signal PDF, we use two-dimensional functions approximated by smooth his- 
tograms obtained from MC simulation. The continuum background PDF is taken as the 
product of shapes in Mb c and AE, which are assumed to be uncorrelated. We use an 
ARGUS [22| parameterization, /(M bc ) oc M^y/l — x 2 exp[— £(1 — x 2 )}, to model the Mb c 
background, with x given by M hc /E heam and £ as a fit parameter. The AE background shape 
is modeled by a normalized second-order polynomial whose coefficients are fit parameters. 

Figure [2] illustrates the fits for the B yields in a baryon-antibaryon mass region below 2.8 
GeV/c 2 , which we refer to as the threshold-mass-enhanced region. The M^ c distributions 
(with -0.135 GeV< AE < 0.074 GeV for tt°/7 modes and \AE\ < 0.03 GeV for the vr" 
mode), and the AE distributions (with M^ c > 5.27 GeV/c 2 ) for the pAj, pAir and pAn~ 
modes are shown. The solid curves show the projections of the fit results. The B yields 
are 98 1]J, 56 jig , and 129 with statistical significances of 14.3, 9.5, and 18.9 standard 
deviations for the pA^, pAir , and pAn~ modes, respectively. The significance is defined as 



21n(L /L max ), where L and L max are the likelihood values returned by the fit with the 
signal yield fixed to zero and at its best fit value. 

Figure [3] shows the differential branching fractions of B + — > pA'-f, B + — > pAir and 
B° — > pAn~ as a function of baryon pair mass, where the branching fractions are obtained 
by correcting the fitted B yields for the mass-dependent efficiencies estimated from MC sim- 
ulation for each mode. Systematic uncertainties are determined using high-statistics control 
data samples. The tracking efficiency is measured with fully and partially reconstructed D* 
samples. For proton identification, we use a A — > pn~ sample, while for K/tt identification 
we use a D* + — > D tt + , D° — ► K~n + sample. The average efficiency difference for particle 
identification (PID) between data and MC has been corrected to obtain the final branching 
fraction measurements. The corrections are about 8%, 8%, and 14% for the pA^, pAn° 
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FIG. 2: The Ai2 and My, c distributions for (a)(b) pA'y, (c)(d) pA-K and (e)(f) pAir~ modes with the 
requirement of baryon-antibaryon mass < 2.8 GeV/c 2 . The solid curve represents the fit projection, 
which is the sum of signal (dash-dotted peak) and background (dashed curve) estimations. 



and pAir~ modes, respectively. The uncertainties associated with the PID corrections are 
estimated to be 2% for protons and 1% for charged pions. For A reconstruction, we have 
an additional uncertainty of 2.5% on the efficiency for tracks displaced from the interaction 
point. This is determined from the difference between A proper time distributions for data 
and MC simulation. There is also a 1.2% error associated with the A mass selection and a 
0.5% error for the A vertex selection. Summing the errors for A reconstruction in quadra- 
ture, we obtain a systematic error of 2.8%. A 2.2% uncertainty for the photon detection 
is determined from radiative Bhabha events. For the 7r° and rj vetoes, we compare the fit 
results with and without the vetoes; the difference in the branching fraction is 0.5%, which 
is taken as the associated systematic error. The uncertainty in 7r° reconstruction is studied 
with D — > Kir and D — > Ktttc samples. The TZ continuum suppression uncertainty is 
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FIG. 3: Differential branching fractions for (a) pAj, (b) pAn and (c) pAir~ modes as a function of 
baryon-antibaryon pair mass. The shaded distribution shows the expectation from a phase-space 
MC simulation. The theoretical predicted curves from Ref. [11] for the pA'y mode and from Ref. 
[10] for the pAn~ mode are overlaid for comparison. The area of the shaded distributions and areas 
under the theoretical curves are scaled to match the measured branching fractions from data. The 
uncertainties are statistical only. 



estimated from B — ► Dtt, D — > Kgir control samples, which have topologically similar final 
states. The determined efficiencies near threshold contribute an error due to the binning ef- 
fect in M p ^. Using the generated MC samples, we vary the bin size to estimate this effect. A 
systematic uncertainty in the fit yield is determined by applying different signal/background 
PDFs and by varying the parameters of the signal and background PDFs by one standard 
deviation. The pAn° mode has a bigger fitting uncertainty due to a larger fluctuation in the 
lower AE side. The error on the number of BB pairs is 1.3%, where we assume that the 
branching fractions of T(4S) to neutral and charged BB pairs are equal. The systematic 
uncertainties for each decay channel are summarized in Table HJ where correlated errors are 
added linearly within each item, and then uncorrelated items are combined in quadrature. 
The total systematic uncertainties are 9.0%, 11.1% and 9.0% for the PA7, pA-zr and pA.n~ 
modes, respectively. 

Table HT1 gives the measured branching fractions for different M p x mass bins. We sum these 



7 



TABLE I: Systematic uncertainties of the branching fraction for each decay channel. 



Source 


pA*y 


pA7T° 


pAn 


Tracking 


4.9 


% 


4.7% 


5.8% 


Proton Identification 


4.0 


% 


4.0 


% 


4.0% 


K/-7T Identification 






- 




1.0% 


BR of A -> pn~ 


0.8 


% 


0.8 


% 


0.8% 


A selection 


2.8 


% 


2.8 


% 


2.8% 


Photon reconstruction 


2.2% 






- 


7T° and 7] veto 


0.5% 








7T° reconstruction 






4.0 


% 




Likelihood Ratio Selection (TV) 


2.5 


% 


4.0 


% 


4.0% 


Modeling and MC statistical error 


3.9 


% 


3.3 


% 


2.0% 


Fitting 


2.2% 


5.6 


% 


1.0% 


Number of BB pairs 


1.3 


% 


1.3 


% 


1.3% 


Total 


9.0 


% 


11.1% 


9.0% 



partial branching fractions to obtain B(B + -> pAj) = (2.45to;38 ± °- 22 ) x 10 ~ 6 > B ( B+ ~> 
pAn ) = (3.00lg:B3 ± 0.33) x 10" 6 , and B(B° -> pAn') = (3.23l£i ± 0.29) x 1CT 6 . These 
values are in good agreement with our previous measurements [a, |8J and supersede them. 
Note that the results include the first observation of B + — > pAir . The ratio of B(B + — > 
pAn°) /B(B° — > pAn~) is 0.93i jg ± 0.09, which is larger than the theoretical prediction of 
0.5. However, one cannot rule out the naive factorization picture with current statistics. The 
shapes of the near threshold peaks can be compared with theoretical predictions [l(J U|, as 
shown in Fig. [3J This comparison is useful for validating (and possibly modifying) theoretical 
models. 



TABLE II: Measured branching fractions 0(10 



for each M ^ bin. 



M pA (GeV/c 2 ) 



pA'y pAir pAn 



threshold - 2.2 
2.2 - 2.4 
2.4-2.6 
2.6 - 2.8 
2.8 - 3.4 
3.4-4.0 
4.0-4.6 



4.6 - M 



pA— lim 



1 02 +0 - 18 
1 - uz -0.16 

Q6+ - 18 
u ' yD -0.16 



0.02 
0.04 
0.03 
0.10 



+0.10 
0.09 

+0.08 
0.08 

+0.13 
0.11 

+0.15 



-0.10 

0.26j£3 

-0.19 



o.oi; 



0.18 



0.21 
-0.17 
0.20 
0.16 
+0.21 
0.17 
0.16 
0.13 
+0.18 
0.18 
0.19 
0.16 
+0.19 
0.16 
n qv+0.33 
U-- 3 ' -0.31 



0.75 
0.54 
0.47 
0.20 
0.14 
0.30 
0.22 



0-88±g : « 
0.86^ 

-0.11 
-0.09 
+0.10 
-0.08 
+0.11 
-0.09 
+0.06 
-0.06 
+0.11 
-0.10 



0.38 
0.22 
0.33 
0.04 
0.23 



below 2.8 
full region 



2.04 
2.45 



0.28 
0.26 
0.44 
-0.38 



1.97 
3.00 



0.39 
0.32 
0.61 
-0.53 



9 o 4 +0.25 
z -°^-0.22 
q 00+0.33 

°-^°-0.29 



We also study the two-body intermediate decays B° — ► pE* , B° — > A°A, B + — > pS*°, 
and B + — > A + A, where the and A 0,+ are reconstructed in the £*- *° —> An~'° and 
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A°' + — > p7i°' + channels, respectively. The selection criteria are 1.30 GeV/c 2 < Mj^ n -,o < 1.45 
GeV/c 2 and M p7T o,+ < 1.40 GeV/c 2 . No significant signals are found in these decay chains. 
We observe 34, 50, 32 and 43 events in the signal region; the expected number of background 
events are 36.9 ± 1.5, 51.8 ± 1.8, 34.0 ± 1.3 and 41.8 ± 1.2 for B° -> p£*~, B° A°A, 
B + — > and B + — > A + A, respectively. We set upper limits on the branching fractions at 
the 90% confidence level using the methods described in Refs. 23|, |24|, where the systematic 
uncertainty is taken into account. The results are summarized in Table UTTl 

In the low mass region below 2.8 GeV/c 2 , we study the proton angular distribution of 
the baryon-antibaryon pair system. The angle 6 P is defined as the angle between the pro- 
ton direction and the meson (photon) direction in the baryon-antibaryon pair rest frame. 
Figure H] shows the differential branching fractions as a function of cos6 p . We define the 
angular asymmetry as Aq = f—^— , where Br + and Br_ stand for the measured branch- 
ing fractions with cos8 p > and cos6 p < 0, respectively. The angular asymmetries are 
determined to be 0.29 ± 0.14(stat.) ± 0.03(syst.), -0.16 ± 0.18(stat.) ± 0.03(syst.), and 
—0.41 ± 0.11(stat.) ± 0.03(syst.) for the pAj, pAn , and pAn~ modes, respectively. A 
systematic error, ~ 0.03, is determined by studying low momentum A reconstruction in 
different angular regions, and by checking the B + — > J/tpK + (J/ip — > sample and 

the continuum background of B + — > ppK + where a null asymmetry is expected. 

Since Aq is not consistent with zero for B° — > pAir~ , the peak near threshold cannot 
be described by a single resonant state [25j | . The opposite slopes in the distributions for 
the pA'y and pAn" modes indicate that the pAj decay agrees well with the short- distance 
b — > sj picture while the pAir~ mode disagrees with the short-distance b —>■ sg description, 
where g stands for a hard gluon. The low mass peaking structure in M p x implies that p and 
A are moving in parallel in the B meson rest frame. One can look for correlations using the 
angular distributions in Fig. HI The negative slope for the pA-K~ mode in Fig. implies 
that the proton moves faster and the A moves slower. In other words, the s quark from b 
decay is not as energetic as expected. Disagreement between data and the short-distance 
description has already been found in the decay B + — > ppK + [8]. One possible explanation 
is the contribution of long-distance effects. 

Another interesting feature of B decays with a A in the final state is the possibility of using 
the A as a helicity analyzer of the s quark in order to check the left-handedness of b — > s 
weak decays. We modify the unbinned likelihood fit in order to simultaneously estimate 
the anisotropy parameter of the secondary proton from A decays. The parameterization 
is 1 + a cos 6, where a is the parameter and 9 is the angle between the secondary proton 
momentum and the direction opposite to the B momentum in the A rest frame. Note that 
the anisotropy parameter a is identical for both A and A. The measured values are — 0.57± 
0.33(statO±0.10(syst.), _-0.27±0.33(stat.)±0.10(syst.), and -0.28±0.21(stat.)±0.10(syst.) 
for the PA7, pAw and pAiT~ modes, respectively. The average A energies in the B rest frame 
are determined to be 1.92 GeV, 1.85 GeV, and 1.78 GeV with standard deviations of 0.33 
GeV, 0.36 GeV, and 0.40 GeV for the pAj, pAir and pAir~ modes, respectively. Figure [5] 
shows the measured anisotropy parameters for different decay modes and compares the 
results with the prediction of the Standard Model 12] as a function of A energy. They are 
consistent within errors. The value of a obtained for the pAir~ mode also agrees well with 



the theoretical prediction in Ref. [10j. The systematic uncertainty in a is included in the 
plot and is about 0.10. This is estimated by varying various selection cuts; the dominant 
effect is the efficiency change near the cos# ~ 1 region, where the detection efficiency for 
slow pions is rapidly changing. 
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We also measure the charge asymmetry as Acp= (Nb — Ni)/(Nf, + iVg) for these modes, 
where b stands for the quark flavor of the B meson. The results are included in Table ITO1 The 
measured charge asymmetries are consistent with zero within their statistical uncertainties. 
The systematic uncertainty is assigned by the measured asymmetry of the background events 
in the candidate region. 
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FIG. 4: Differential branching fractions vs. cos6 p for (a) pA'y, (b) pAir and and (c) pAir~ modes in 
the region near threshold (baryon-antibaryon mass < 2.8 GeV/c 2 ). The uncertainties are statistical 
only. 

In summary, using 449 xlO 6 BB events, we measure the mass and angular distributions of 
the baryon-antibaryon pair system near threshold for the pA^, pAir and pAir~ baryonic B 
decay modes. We report the observation of B + — > pAn° with a branching fraction (3.00jTQ;| 3 ± 
0.33) x 10~ 6 and a low pA mass peak near threshold. The measured branching fractions for 
B + — > pA'-f and B° — > pAn^ are in good agreement with our previous measurements @]. 
The different proton polar angular distributions for the pA'j and pA^~ modes indicate a 
difference between b — > S7 and b — >• sg decays. The anisotropy parameters a from A decays 
agree with theoretical predictions within errors. We also search for intermediate two-body 
decays and find no significant signals. We set upper limits on their branching fractions at 



the 90% confidence level. Some suppression factors [26| for the charmless baryonic two-body 



decays should be considered under the present theoretical framework and understanding the 
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TABLE III: Summary of the measured results for B + — > pA^, pAir® and B° — > j>A"7r~. Y is the 
fitted signal or upper limit at 90% confidence, a is the statistical significance, B is the branching 
fraction, Aq is the angular asymmetry and Aqp is the charge asymmetry. 



Mode 


Y a B (10~ 6 ) 


A e 


A C p 


B + — »■ _pA7 


114+{^ 14.5 2.45±g;|| ± 0.22 


0.29 ±0.14 ±0.03 


0.17 ±0.16 ±0.05 


5+ pAvr 


8911? 10 - 2 3.00lg;g ± 0.33 


-0.16 ±0.18 ±0.03 


0.01 ±0.17 ±0.04 


5+ -► 


< 11.3 - < 0.47 






5+ -► A+A 


< 15.9 - < 0.82 






5° pAvr" 


1781^ 20 -° 3 - 23 ±a! ± °- 29 


-0.41 ±0.11 ±0.03 


-0.02 ±0.10 ±0.03 


B° ^pt*- 


< 10.9 - < 0.26 






B° A°A 


< 15.9 - <0.93 







mechanism of the threshold enhancement might be the key to determine the two-body decay 
rates . 
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